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A  control  system  was  designed  to  attenuate  vertical 
accelerations  for  the  XR-3  captured  air  bubble  type  surface 
effect  ship  using  linear  regulator  techniques  applied  to  the 
simplified  nonlinear  equations  of  motion.   A  pressure  lift- 
only  model  was  used  to  represent  the  craft  vertical  heave 
motion  and  was  linearized  around  the  steady  state  operating 
point.   Model  validation  was  obtained  through  analysis  of 
the  frequency  spectrum.   State  variable  feedback  was   used  to 
determine  a  set  of  optimal  control  gains  that  would  reduce 
the  magnitude  of  the  heave  acceleration  during  operation 
under  simulated  sea  input  conditions. 
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I,   INTRODUCTION 

An  important  design  consideration  for  the  high  speed 
captured  air  bubble  type  surface  effect  ship  is  the  vertical 
acceleration  associated  with  operation  in  a  random  ocean 
environment  [Ref,  1],   The  vertical  motion,  heave,  is  caused 
by  waves  changing  the  volume  of  the  plenum  producing  pres- 
sure variations  that  result  in  vertical  craft  accelerations. 

Studies  of  conventional  hull  ships  show  that  the  strongest 
vertical  plane  accelerations  are  experienced  at  encounter 
frequencies  close  to  the  ship's  natural  frequency.   Reference 
2  states  that  "The  most  severe  and  objectionable  motions  (wet 
decks,  slamming,  high  accelerations)  result  from  the  responses 
to  the  components  that  are  near  synchronism  with  the  ship's 
natural  pitching  and  heaving  periods".   Encounter  frequency 
is  a  function  of  wave  frequency,  craft  speed  and  direction. 
Since  most  of  the  weight  of  the  surface  effect  ship  is  sup- 
ported by  a  bubble  of  air,  large  accelerations  can  occur  that 
affect  habitability  and  sea  keeping  ability  while  operating 
in  high  seas . 

Methods  of  attenuating  vertical  motion  have  been  the 
subject  of  past  studies.   A  conventional  design  by  Rohr  Marine 
Inc.  [Ref.  3]  utilized  accelerometer  output  to  monitor  verti- 
cal motion  variations.   These  signals  are  passed  through 
compensation  circuits,  filtered  and  applied  to  actuators 


which  control  air  flow  rates  into  the  plenum.   State  variable 
techniques  are  applied  to  describe  the  system  and  control 
functions.   The  physical  constraints  of  the  actuators  are 
considered  with  regard  to  their  inherent  nonlinearities  and 
limitations,   The  sea  input  is  modeled  using  a  wave  height 
power  spectral  density  procedure. 

Through  the  use  of  a  towing  tank  model  and  data  gathered 
from  the  100B  test  craft,  an  investigation  was  conducted  to 
measure  the  characteristic  response  of  the  SES  while  operating 
in  an  actual  sea  state  environment  [Ref.  4],   The  model  was 
subjected  to  simulated  sea  state  inputs  of  various  severity  in 
an  effort  to  establish  the  linear  characteristics  of  the 
vertical  motion  and  determine  the  damped  natural  frequency  of 
response.   A  series  of  curves  showing  response  amplitude  for 
both  C.G.  accelerations  and  pitch  angle  versus  encounter  fre- 
quency were  developed  for  a  range  of  speed  and  sea  state 
operation . 

In  Ref.  1,  the  frequency  response  approach  was  used  in  an 
analysis  which  examined  the  effects  of  changing  plenum  air 
flow  rates  on  center  of  gravity  accelerations  using  the  six 
degree  of  freedom  simulation  program  of  the  XR-3  [Ref.  5]. 
This  study  developed  a  linear  simplified  model  in  close  agree- 
ment with  the  6  D.O.F.  model  and  in  addition  used  signifi- 
cantly less  computer  time, 

In  Ref.  6,  Boggio  conducted  a  study  on  the  effect  of  the 
use  of  a  flexible  membrane  in  the  plenum  chamber  itself. 
Using  the  loads  and  motions  program  for  the  XR~3 ,  the  study 


reduced  computation  time  through  the  inclusion  of  pressure 
rate  and  volume  rate  equations.   Membrane  insertion  served  to 
effectively  reduce  the  magnitude  of  pressure  variations  in 
the  plenum  thereby  attenuating  vertical  motions. 

Grant  [Ref,  7]  used  simple  models  for  the  control  devices 
applied  to  the  loads  and  motions  program  of  the  100B  SES 
craft.   From  monitored  plenum  pressure  variations  and  heave 
accelerations,  control  signals  were  developed  by  filter  de- 
sign methods ,  to  vary  the  fan  speed  and  vent  louver  openings 
in  order  to  reduce  the  vertical  accelerations  and  pitch  motion. 

This  thesis  applies  the  linear  regulator  design  technique 
to  determine  a  set  of  optimal  state  feedback  gains  to  control 
the  air  flow  rate.   The  equations  of  motion  used  in  the  design 
were  derived  from  a  simplified  pressure  lift  only  model  lin- 
earized around  the  steady  state  operating  point.   Design  re- 
sults were  evaluated  using  both  the  linear  and  nonlinear 
models  . 


II.   PRESENTATION  OF  THE  SYSTEM 

Two  major  forces  which  support  the  weight  of  the  surface 
effect  ship  are  the  buoyancy  and  plenum  pressure  forces. 
The  buoyancy  forces  (F,)  are  exerted  on  the  craft  by  dis- 
placement of  water  by  the  underwater  volume  of  the  rigid 
sidewalls.   The  plenum  pressure  lift  force  (Fp)  is  exerted 
on  the  craft  by  the  captured  cushion  of  air.   Flexible  bow 
and  stern  seals  along  with  the  rigid  sidewalls  and  lift  fans 
form  the  captured  air  bubble.   Relative  contributions  to 
craft  vertical  plane  motion  by  FT  and  Fp  depend  upon  the 
encounter  frequency  (oj  ).   At  low  frequencies,  FT  responds 

to  the  rise  and  fall  of  the  sea.   At  high  oo  ,  FT  attenuates 

°        e   L 

because  incremental  volume  changes  of  the  immersed  hull  are 
reduced.   The  Fp  force  responds  inversly  to  changes  in  plenum 
volume  (V_)   according  to  the  adiabatic  law  of  pressure 

D 

variations  within  the  plenum  chamber  [Ref.  1].   Volumetric 
changes  in  the  plenum  depend  upon  the  sea  state  and  encounter 
frequency.   As  is  shown  in  Ref.  1,  the  effect  of  the  volume 
changes  produces  an  increase  in  the  Fp  force  as  the  encounter 
frequency  is  increased. 

In  the  system  developed  in  this  study,  the  buoyancy  terms 
have  been  omitted  in  an  effort  to  simplify  the  system  further, 
This  simplified  model  is  referred  to  as  the  "pressure  lift 
only"  system.   Transfer  functions  of  acceleration,  volumetric 
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plenum  air  mass  to  input  waveheight  and  system  state  equa- 
tions were  developed  from  the  linearized  form  of  the  pres- 
sure lift  only  system. 

The  basic  assumptions  which  form  the  basis  for  the  pres- 
sure lift  only  model  are  presented  below: 

1.  The  center  of  pressure  (C  )  is  directly  under  the 

P 

center  of  gravity  (C.G,), 

2.  The  sidewalls  are  of  uniform  cross  section  and  are 
symmetrical  about  the  C.G. 

3.  The  effect  of  pitch  moments  on  the  vertical  plane 
forces  are  neglected. 

4-,   The  only  variable  lift  forces  are  those  forces 
caused  by  changes  in  plenum  pressure  (P,  ). 

5.  The  craft  is  at  a  constant  speed. 

6.  The  air  leakage  rate  from  the  rear  seal  is  constant. 

7 .  The  craft  weight  is  supported  by  two  component 
forces  in  the  following  approximate  proportions ; 

Sidewalls:   20%  (fixed  force) 

Captured  Air  Bubble:   8  0%  (variable  force) 

8.  Encounter  frequencies  of  interest  lie  within  the 
range  of  17<oj  <100  radians/sec. 

A,   EQUATIONS  OF  MOTION 

The  development  of  the  nonlinear  simplified  model  for 
vertical  motion  is  presented  in  Ref .  1  and  the  block  diagram 
appears  in  Appendix  A.   The  development  of  the  pressure  lift 
only  nonlinear  model  follows  in  a  similar  manner.   The  equa- 
tions of  motion  differ  in  that  the  F.  force  is  considered 

constant . 
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The  application  of  the  given  assumptions  to  the  Newtonian 
Laws  for  vertical  motion  of  a  rigid  body  system  leads  to  the 
following  equations; 

1)   Orifice  Leakage  Rate 


=  C  A1      t/  — 

>ut    n  1   W  p 


C  A.,  \l — Z  ££~  (2.A.1) 

^•out    n  1   W  p        sec 


2)   Fan  Map  Input  Rate 


ft3 
q.   =  nCQ.(o)  -  K  P.]    jj-  (2. A. 2) 

^in       i        q  b     sec 


3)   Absolute  Plenum  Pressure 


M,    y     1K 
P.  =  P   (  tt^-  )      iSI  (2. A. 3) 

b    a    Vbpa       ft2 


4)   Plenum  Air  Flow  Rate 


&  =  p   (q.   -  q   O    —  (2. A. 4) 

b    a   Hm   Hout     sec 


5)   Heave  Acceleration 


Z  "  (M  ""  IT  -  IT5   iic  U.A.b; 


6)   Heave  Velocity 

Z  =    f'i   dt  —  (2. A. 6) 

./  sec  s 
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7)   Heave  Displacement 


Z  =  f     Z  dt         ft  (2, A. 7) 


8)   Plenum  Pressure  Lift  Force 


Fp  =  AbPb  lbs  (2. A. 8) 


9)   Buoyancy  Lift  Force 


Fl  =  C2As1dpg]        lbs  (2. A. 9) 


where  1 ,  is  steady  state  draft 


10)   Plenum  Gage  Pressure 


PK  =  (P.  -  P  )        i£|  (2. A. 10) 

b      b     a        ft2 


11)   Draft 


1 ,  =  Z  +  Z  ft  (2. A. 11) 

d        s 


The  system  parameters  and  constants  are  listed  as  follows : 


1)   Adiabatic  Process  Coefficient 


Y  =  1.4 
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2)   Leakage  Area 


A±    =  YL  ZL  =  ,438     ft2 


3)   Leakage  Orifice  Coefficient 


C   =  .90 
n 


4)   Air  Density 


p   =  .002378 
ya 


5)   Atmospheric  Pressure 


P   =  2116 
a 


6)   Plenum  Area 


A^  =  200  ft2 

D 


7)   Empty  Plenum  Volume 


V   =  383  ft3 

n 


8)   Craft  Weight 


W  =  6007  lbs 
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9)   Area  of  Each  Immersed  Sidewall 


A  =  75/4  =  18.75      ft2 


10)   Draft  Initial  Condition 


Ld(o)  =  .5  ft 


11)   Density  of  Water 


p  =  1.99 


12)   Gravity  Acceleration  Constant 


G  =  32.174  ft/sec2 


13)   Vertical  Center  of  Gravity  Location  Above  Keel 
Line 


Z   =  2.5  ft 

s 


140   Number  of  Fans 


n  =  5 


15)   Steady  State  Fan  Output 


3 
Q.  (o)  =  35  ft  /sec 
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16)       Fan   Map    Slope 


K      =    1.0 

q 
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B.   LINEAR  MODEL 

Prior  to  defining  the  system,  states ,  the  non  linear 
equations  were  linearized  around  their  steady  state  operating 
conditions  by  application  of  a  Taylor  series  expansion.   The 
details  of  the  linearization  process  is  shown  in  Appendix  A. 

The  resulting  simplified  pressure  only  signal  flow  graph 
for  the  linearized  system  is  shown  in  Figure  2 . 

1 .   Transfer  Function  Development 

The  system  graphical  determinent  is; 


A  =  1  -  (L   +  L2)  +  h1   L2  (2.B.1) 


A  =  1  -  (Lx  +  L2)  (2.3.2) 


where 


L   =  s"1s"1(-k8)(-k2)(k3)(-k1+) 


=  -s'VsVs  (2-B-3) 


L0  =  s"1(k7)(kJ(-k_) 
2         7    a    6 


=  -3-1k3kgk7  (2.B.U) 


The  transmission  gain  from  the  control  input  to  the  accel- 
eration output  for  calm  water  operation  is ; 
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»    P1A1 

3z    =    -A^      \  (2.B.5) 


where 


Pl    =    ^k3k4k7k10S"1  (2.B.6) 


and 


Ax    =    1  (2.B.7) 


and 


aT7    =    0  (2.B.8) 

W 


Substituting  equations  (2.B,2),  (2.B.6),  (2.B.7)  into 
equation  (2.B.5)  yields; 


—  k-k.  Kr7J<C,  n  S  3 

3"  =    3  ±    I    10 SL. (2.3.9) 

l-[-s"zk2k3k4k8)+(-s"xk3k6k7)] 


Reducing  further; 


«T  =   23  "  ?  1Q *  (2.B.10) 

s  +k3k6ky  s+k2k3k4k8 


The  transmission  gain  from  the  sinusoidal  sea  generator  input 
to  the  acceleration  output  with  no  control  signal  is; 

P2A 

3"  =  — —   3TJ  (2.B.11) 

z     A     W 
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where : 


P2  =  "klk2k3k4  (2,8.12) 


and 


Ax  =  1  (2.B.13) 


and 


3U  =  0  (2.B.14) 


Appropriate  substitution  into  equation  (2.B.11) 


yields 


3'  -      2  3  4 


1-C(s"  k2k3kuk8)+(-s"  k3kgk7)] 


Reducing  further,  the  transfer  function  is 

I  =  -;ik2k3^  s2 (2.B.15) 

s  +k3k6k7s+k2k3kuk8 

For  the  condition  that  control  input  is  applied,  the 
acceleration  output  has  the  following  transform? 
-k„k,.k„k.  „s3 


-3  4JV7  N10   u 


Z    s  +k3k6k7S+k2k3k4k8 

k1k0k„k„s2    9W 

„   l±    I    3    4 , C2.B.16) 

s  +k3k6k7s+k2k3k4k8 
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Following  the  same  procedure  for  the  applied  control 
input,  the  volumetric  plenum  air  mass  output  has  the 
following  transfer  function: 

A3  =  1  :<*   L±  (2,B,17) 


P3  =  k1Qs"1  (2.B.18) 


A2  =  1  (2.B.19) 


P^  =  -k1k2k3k6  (2.B.20) 

3Vm  =  y^^^8""3"1    3u  (2.B.21) 

s  «*k «x nKn s  ""k 0K „ k,, k_ 
3  5/    2348 


Equation  (2.B.15)  is  of  the  form; 


OUTPUT      -as2 


INPUT  "  ;2..   . 

S  +  DS  +  C 


and  with  no  control  action,  becomes 

2 


y(s)  as' 

S     +DS+C 


m  _  d 

WTsT    "  7T7~ 

3     +DS+C 


1 — >  <^ 


(2.B.22) 


(2.B.23) 


where 

y(s)    =    3Z  (2,B.2"4) 

and 


(2.B.25) 
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where  the  coefficients  consist  of  the  following  parameters ) 


a-=-k1k2k3k4  (2.B.26) 


b  =  k0kKk7  (2.B.27) 

3  6  7 


c  =  k2k3k4k8  (2.B.28) 


d  =.-k,k0k0kc  (2.B.29) 

12  3  6 


The  values  of  the  coefficients,  k  ,  for  the  XR-3  craft  are 
listed  in  Table  I  and  were  obtained  from  the  results  of  the 
linearization  process  of  Appendix  A, 

Special  attention  should  be  given  to  the  coefficient 
that  contains  k  which  is  a  sine  function  of  the  form: 

2W  V         co.L 

k   =  [(— E-£)  Sin  (*&=&>]  (2.B.30) 

1       co.  I  vw 

l 

k  may  be  simplified  by  an  approximation  detailed  in  Appendix 
A  of  Ref.  3.   The  simplification  is  based  upon  the  fact  that 
for  ahead  seas,  the  incident  frequency,  ok  is-, 

/     Vs 

1-/  l+i+C-JOw 

^  e  (2.B.31) 

r 

■ 

For  convenience  in  this  study,  the  value  of  the  ships  velocity, 

V    is  set  equal  to  the  gravitational  constant,  G.   This  value 
s  ' 

of  V  ,  when  converted  to  knots,  is  19.3  knots  which  is  within 
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TABLE  I 


SYSTEM  COEFFICIENT  (k  )  VALUES 

n 


n  FORMULA 


XR-3    VALUE 


(2\\     .     ALP. 
i  w 


"2 

V- CO) 

k3 

YPb(O) 

\ 

Ab 

M 

nk      + 

q 

C   A, 
n   1    j 

k6 

/    '\ 

P 
a 

2PT(0) 
D 

ky 

1 

V    (0) 

m 

k8 

Ab 

.     ,62.8319 


sjLnC  ■  ■;•■■  --0 


2267,782 


3,5336x10 


2996.0 


1,0711 


6,14422 


m, . 


/6  2.  8"319x 

c  ■  l;;  '  ^ 

w 


-3 


3.5052x10 


-3 


200.0 


k    Q         CONSTANT 


1.0 
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the    operating   range   of    the    XR-3,       Continuing; 

1-    /l+i+w 
w      =    — ^ .  .  ■  .  f.  , — I  C2.B.32) 

1  ~2 

It  follows  that 


oo.2  =  go  (2.B.33) 

i     e 


This  approximation  is  valid  for  all  oj  >2.5  radians/sec.   From 

e 


the  definition  of  the  incident  wave, 


cd.2  =  ^  (2.B.34) 


k,  becomes 


'aiLP 

Sin^  2V 


kn  =  A,  — r-^ (2.B.35) 

(  3=  p) 

^  2V 


OJ 


Where  the  horizontal  area  of  the  plenum  chamber  is 


Ab  =  WPLP   '  (2.B.36) 


For  a  speed  of  19.3  knots  and  since 


LT7  =  ^&  C2.B.37) 

e 
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k   simplifies  toj 

•itL 


kl  =  Ab   ;.^lW—  (2.8,38) 


For  the  XR-3 , 

A,  =  200  ft2 

D 


L   =  20  ft 
P 


Th6n  C4  ,62.8319, 

LW 

a  =  k   =  2267,782  — — . — . , (2.B.39) 

1  ,€2.8319, 

From  Equation  (2.B,38),  it  can  be  seen  that  a  null  value  for 

k,  will  occur  for  nLT7  =  L   and  a  peak  value  for  k.,  will  occur 
1  W    p        r  1 

for  nLT7  =  L  /2  where  n  is  an  integer. 

W     p  ° 

2  ,   State  Equation  Development 

The  state  equations  are  of  the  form 

X  =  AX  +  B3u  +  G3W  (2.B.40) 


where  the  state  variable  vector  is  chosen  as 

X  =  (3  ,  3"  ,  SV  )T  (2.B.41) 

—     z '   z '   m 


and  the  output  variable 


y=3=CX  +  F3W  (2.B.42) 


z    p 
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The  states  are  defined  as 


XI  =  3Z  =  (Z~Z   K  Incremental  Heave(Draft)    (2.B.43) 
ss  ' 

X2  =  3Z  =  (Z-Z   ),  Incremental  Heave  Rate      C2.B.44) 
ss  * 

X3  =  3V   =  (V  -V    ) ,  Incremental  Air  Mass 

m     m  mss    yolume  (2.B.45) 


where  V   =  r-2.  (2.B.46) 

m   p3 
a 


with  control  defined  as  a  scaler,  rate  of  incremental  change 
of  air  mass  volume, 


3   =  3V  (2.B.47) 

u     m 
c 


and  where  the  sea  input  disturbance  is  3W. 

From  the  signal  flow  graph  of  Fig.  2, 


3V      =    (3    k,n=kc3P^)s"1  (2.B.48) 

m  u   10      6      b 


X.    =    9      =    s"1^  C2.B.49) 


Xn    =     3*    =    s'1^'  C2.B.50) 

1  z 


X2    =    3Z    =    -    3Pbku  C2.B.51) 
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where 

9Pb=(3Vmk7~k2["^Wkl"k8BZ])k3  C2.B.52) 

Rearranging 

2?b~  3Vmk3k7+klk2k3  9w+k2k3k8  3Z  ( 2 . B . 53 ) 

and   grouping    in  order  of   states 


3P,     =    k.k,kQaZ+kQk„3V    +knk0k03W  (2.B.54) 

b  238  37      m      123 


The  output  term,  y  is 


y  =  3  (2.B.55) 

z 


Substituting  equation  (2,B,5)  into  equation  (2.B.55)  yields, 

y=-kok0k.,kQ3Z-k0kllk„3V  -k.kok0k,  3W  (2.B.56) 

Z  o  4  o      3  4/   m   ±  z  d  4 

Finally,  the  incremental  change  in  the  volumetric  air  mass 
is  written  as 


X3  =  3V  =(3uv-kc3P,  )  (2,Bt57) 

m     10   6   b. 


=    ~kok0kcko3Z~k0kck73V 
2368  367    Tm 


-  k1k2k3k63W+k103u  (2,B.58) 


28 


Equations  (2.B.50),  (2,B,56)  and  C2,B.57)  are  the  states  of 
interest  and  are  grouped  below, 


XI  =  X2 


(2.B.59) 


X2    =^k2k3kl+k8^l-k3k4k7:X3-k1k2k3k43W 


(2.B.60) 


X3  =  -kok0kckQXl-k0kck7X3~k,kok0k„ SW+kn ndu  (2.B.61) 

2  3  6  o     36/     1236     10 


with  the  output  acceleration; 


y    =    -k2k3kl+k8Xl-k2kuk7X3-k1k2k3kl+ 3W 


(2.B.62) 


Thus  the  linear  system  matrices  are; 


A  = 


~k2k3k6k! 


""k3k4k7 
"k3k6k7 


(2.B.63) 


B  = 


(2.B.64) 


r  =  —  K..  K„K^iC|^ 


(2.B.65) 


G  s 


"klk2k3k4 


12  3  6 


(2,B.6B) 
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Computer  simulation  of  the  linear  system  was  accom- 
plished using  the  Naval  Postgraduate  School  program  "IODE  ; 
operated  from  the  CPCMS  terminal  of  the  IBM-360.   See  Appendix 
C  for  IODE  listing.   The  program  ocean  generator  for  dVl   inputs 
a  single  frequency  at  a  time  using  the  special  function  fea- 
ture,  Wavelength  (L.,)  will  change  with  encounter  frequency 

w 

( cu  )  as  shown  in  equation  (2.B.35),   It  should  be  remembered 
that  the  craft  speed  is  constant  at  19,3  knots  for  the  data 
base  of  this  design, 

3 .   The  Characteristic  Equation 

The  characteristic  equation  under  no  control  condi- 
tions of  operation,  is  formed  from  the  A  matrix  (see  Appendix 
A), 


DET(SI-A) 


(2.B.67) 


DET 


2267.782 
13002.8 


-1 

s 

0 


0 

11.248 

64.493s 


(2.B.68) 


Then  the  characteristic  equation  is; 


s3+64,493s2+2267,782s+,5701 


(2.B.69) 


which  factors  into 


-4 


(3  +  2,5414x10    ^)Cs+3  2.2  5+J3  5.04)Cs  +  3  2.2  5-j3  5.04) 


C2.B.70) 
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and  is  shown  in  the  pole^zero  map  of  Fig.  3.   Note  that  the 
dominant  pole  is  a  real  pole  close  to  the  origin,   Also  note 
that  the  complex  pair  will  affect  early  time  response  since 
they  are  relatively  far  removed  from  the  origin. 

A  root  locus  plot  with  no  control  input  is  shown 
in  Fig.  4.   The  plot  covers  a  range  of 

5  <  k6  <  26  (2.B.71) 

Table  II  lists  the  characteristic  equation  roots  for  the 
range 


5  <  kc  <  50  (2.B.72) 

o 


As  can  be  seen  from  equation  (2.B.63),  kfi  is  present 
in  both  the  a-,  and  a  _  system  coefficients.   In  order  to 
perform  the  root  locus  study,  the  characteristic  equation 
was  rearranged  to  show  the  adjustable  parameter  k  *, 

s3+k3k7s2+k2k3kl+k8s+k6(k3k4k7k2]<:3k8  +  k2k3kl+k8k3k7)+k6s2 

(2.B.73a) 
with  parameter  values 


s3+(k  +10.5016)s2t2  2  67,78  2s+.09  28kc  (2,B.73b) 

6  -  b 


From  the  table  and  root  locus  plot,  it  can  be  seen  that  as  kg 
becomes  large ,  the  complex  root  pair  become  real  with  one 
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root  moving  toward  the  origin  and  the  other  migrating  to 
negative  infinity.   Additionally,  the  pole  near  the  origin 
moves  very  little  and  remains  real  throughout  the  excursion 
of  kg.   It  was  observed  that  this  was  very  similar  to  the 
root  locus  plot  in  the  air  flow  study  conducted  in  Ref,  1. 
However,  the  removal  of  the  buoyancy  term  in  this  study 
caused  the  pole  near  the  origin  to  stay  real  whereas  the 
same  pole  in  Ref.  1  became  immaginary  for  large  values  of  ke  . 
The  root  locus  demonstrated  the  high  sensitivity  of  the  com- 
plex pair  to  changes  in  air  flow  rate ,  while  the  pole  near 
the  origin  exhibited  a  very  low  sensitivity  to  changes  in 

k_ .   This  observation  further  demonstrated  that  the  vertical 
6 

oscillatory  characteristics  of  the  ship  depend  upon  changes 

in  k_  . 
o 

As  a  further  check  on  the  effect  of  air  flow  on 

C.G.  acceleration,  an  investigation  of  the  effect  of  changing 

kR  on  the  frequency  spectrum  of  C,G.  accleration  was  conducted 

and  the  results  are  shown  in  Appendix  D,   From  these  results 

it  seemed  appropriate  to  seek  a  control  law  that  would  adjust 

k_  ,  the  air  flow  rate,   In  an  effort  to  gain  a  better  under- 
b 

standing  of  the  effect  of  various  air  flow  rates  on  the  sys- 
tem, a  sensitivity  and  stability  analysis  was  performed  in 
section  IVB,   The  R00TL0  program  used  in  the  above  study  is 
shown  in  Appendix  C, 
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Y  ( -32 , 25  + j  35.04) 


* 


2. 5x10 


X (-32.25-j35.04) 


Figure  3 

POLE/ZERO  PLOT  FOR  THE  CHARACTERISTIC  EQUATION-  WITH  NOMINAL 
AIR  FLOW  RATE  s  +64  .  493s  +2267  .  7 82s+ . 57 01 
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Ill,   VERIFICATION  OF  THE  MODEL 

A,   LINEAR  AND  NONLINEAR  SYSTEM  VALIDATION 

The  validation  of  the  simplified  model  of  a  single  vari- 
ble  pressure  lift  force  was  conducted  by  comparison  of  its 
frequency  response  characteristics  with  those  of  the  model 
presented  in  Ref.  1,   Several  techniques  were  used  to 
accomplish  this, 

1.  Single  frequency  constant  amplitude  sinusoids  were 
applied  to  both  the  nonlinear  system  of  Fig.  1  and 
also  the  linear  system  of  Fig.  2.   The  acceleration 
output  signal  for  the  nonlinear  system  is  shown  in 
Fig.  5.   Its  sinusoidal  form,  in  steady  state  indi- 
cates the  linearity  of  the  system,  which  was  the  case 
throughout  the  frequency  range , 

2.  The  linear  transfer  function  of  center  of  gravity 
acceleration  (Eq.  2.B.15)  was  subjected  to  a  constant 
wave  height  input  using  a  spectrum  of  discrete  en- 
counter frequencies.   The  heave  accelerations  were 
converted  to  dB  and  plotted  in  a  bode  diagram  shown 
in  Fig.  6 

3.  Discrete  fourier  transform  technique,  developed  in 
an  independent  study  by  S,  Carpenter,  [Ref,  9],  was 
also  used  to  test  the  response  of  the  linear  system 
(Fig,  9)  to  an  irregular  sea  input  shown  in  Fig.  7 
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consisting  of  twelve  components  of  frequency  and 
amplitude  from  the  Pierson-Moskowitz  energy  density 
spectrum  given  in  Fig.  8  [Ref.  8,  9]. 
The  results  of  procedures  2  and  3  above  were  compared 
with  the  linear  system  response  of  Ref,  1,   The  magnitudes 
for  the  center  of  gravity  acceleration  showed  approximately 
20%  difference  between  the  pressure  lift  only,  no  control, 
response  curve  (Fig.  6)  and  the  linear  system  of  Ref.  1 
(Fig.  10),   The  general  shape  of  the  curves  are  identical. 
Therefore,  the  pressure  lift  only  model  was  assumed  to  be 
valid  and  useful  for  linear  regulator  control  design.   The 
observed  peaks  and  nulls  of  Figs.  6,  10  are  a  result  of  the 
sine  function  numerator  coefficient,  k  ,  as  described  in 
Section  III.   The  relationship  between  plenum  length  (L  ) 
and  wavelength  (L  )  expressed  in  k, ,  determine  where,  in 

W  J. 

the  encounter  frequency  spectrum,  the  nulls  and  peaks  occur 
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SYSTEM  STABILITY  TESTS 
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IV,   THE  CONTROL  LAW 

A.   LINEAR  REGULATOR  DESIGN  APPROACH 

The  linear  pressure  lift  only  model  is  described  in 
state  equation  form; 
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with  the  cost  function  of  the  form 


J  =  1/2 
0 


/ 


(XTQX  +  8UTR3U) 


(4. A. 2) 


where  the  final  value  error  is  not  included.   The  system  may 
be  designed  using  the  linear  regulator  technique  with  the 
given  quadratic  cost  function.   An  optimal  control  was 
found  using  both  the  transient  and  steady  state  solutions 
of  the  Riccati  equation.   The  control  law 


*     -1  T 

3U   =?R   B  KX 

-K     ;>-   ft  — 


(4. A. 3) 


where  the  K  matrix  was  found  from  the  solution  of  the  Riccati 
equation 
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K  *  -KA  -  ATK  -  Q  +  K  B  K^BTK 
~     "•<     *i     ft  n      n 


(4.A-4) 


with  the  appropriate  choice  of  the  Q  and  R  coefficients,  the 
feedback  gains  for  the  system  were  calculated.   Unfortunately, 
the  first  choice  for  Q  and  R  may  not  yield  control  gains  which 
are  physically  realizable  or  practical.   For  this  reason,  it 
is  appropriate  to  examine  the  effects  of  the  control  gains  on 
the  pole  location  in  the  s  plane  for  the  controlled  system. 
Such  a  study  was  done  and  is  presented  in  part  B  of  this 
section.   The  choices  for  Q  and  R  were  derived  by  trial  and 
error.   However,  the  design  criterion  was  to  reduce  this  mag- 
nitude of  the  frequency  response  curves  through  the  study  of 
various  Bode  plots.   To  this  end,  many  Q  and  R  values  were 
tried . 

The  Riccati  equation  computations  were  done  using  a  com- 
puter program  documented  in  Ref,  10  and  modified  by  B. 
DesJardins  in  Ref.  11.   Appendix  C  illustrates  the  input/ 
output  format  of  this  program.   Based  upon  the  above  design 
criterion,  the  final  weighting  matrix  chosen  was; 


Q  = 


10  0  0 
0  103  0 
0       0       .5 


(4. A, 5) 


and 


R  =  .03 
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The  optimal  control  gains  obtained  using  the  above  Q  and  R 
matrixes  were; 

gll  =  -1.689xl03 

g12  =  l,19xl02  (4. A. 6) 

g13  =  -1.83X101 

These  gains  were  multiplied  by  their  respective  states  and 
applied  to  the  state  equations  (2.B.59,  61,  62).   Under  state 
feedback  control,  the  resultant  output  acceleration  frequency 
magnitude  was  reduced  by  13%  from  the  peak  no  control  magni- 
tude.  Table  III  lists  the  peak  lobe  values  of  the  C.G, 
Acceleration/wave  height  magnitudes  versus  frequency  for  the 
nonlinear  and  linear  systems.   Each  system  was  exercised 
under  control  and  no  control  conditions.   The  steady  state 
sinusoidal  output  for  the  linear  model,  under  control  and 
no  control  conditions  is  shown  in  Table  IV    for  ^e  =  t+5  rad/sec, 
where 

3u  =  SllX1  +  §12X2  +  Sl3X3  (4-A-7) 

and 

y  =  3Z  (4. A. 8) 

The  control  signal  for  the  nonlinear  system  was  generated  by 

using  the  linear  control  low  coefficient  found  using  linear 

regulator  theory.   The  reduction  in  acceleration  was  not 

uniform  throughout  the  frequency  range.   At  law  frequencies, 

to  <17  rad/sec,  control  gains  caused  an  increase  in  C.G.  acceleration. 
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As  «e  approached  100  rad/sec3  C.G,  acceleration  approached  but 
never  exceeded  the  no  control  condition.   Effective  control 
was  realized  in  the  frequency  range  17 <W  <100  radians/sec. 


e 


B,   SENSITIVITY  AND  STABILITY 

Since  a  set  of  optimal  gains  were  found  for  the  system, 
it  remained  to  test  the  system  under  feedback  control  and 
observe  the  sensitivity  and  stability  for  acceptability. 

The  locus  of  roots  of  the  characteristic  equation  for 
the  pressure  lift  only  system  without  state  feedback  control 
was  shown  in  Section  III,   The  objective  of  state  feedback 
control  is  to  adjust  the  position  of  the  complex  pair  of 
roots  in  such  a  way  as  to  reduce  center  of  gravity  accelera- 
tions.  For  the  open  control  loop,  the  transfer  function 
was  shown  (Eq.  2,B,22)  to  be 


2 
dZ      _     -as'  ,,  o   . 

w    -      rr  7  (4.B.D 

s  +bs+c 


Addition  of  state  feedback  will  change  the  value  of  b  and  c. 
Note  that  if  the  k,  sj;nc  function  null  effects  are  removed 
from  the  numerator,  thejia  =  c.   It  follows  that  the  magnitude 
of  the  Bode  plot  will  reach  unity  value  at  high  frequency 
regardless  of  changes  in  b  and  c,   In  other  words,  the  sys- 
tem has  high  pass  frequency  characteristics.   Thus,  the 
results  of  increasing  the  value  of  c  will  be  a  reduction  of 
the  magnitude  plot  in  the  lower  range  of  frequencies  of 
interest.   Also,  note  that  changes  in  the  b  coefficient  will 
affect  the  resonant  peaking  at  the  corner  frequency, 
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TABLE  III 

PEAK  LOBE  VALUES  OP  FREQUENCY  RESPONSE,  PRESSURE  LIST  ONLY 
LINEAR  S  NONLINEAR  MODELS  UNDER  CONTROL  AND  NO  CONTROL 

CONDITIONS 


CO 

e 

NONLINEAR 

LINEAR 

NO  CONTROL 

CONTROL 

NO  CONTROL 

CONTROL 

6 

18,6 

26,5 

18,5 

27,8 

7 

18,5 

25.2 

18.6 

26,6 

8 

15.7 

20,3 

15.7 

21.4 

9 

9.7 

12,7 

9,8 

12.7 

14 

42  .2 

42.2 

42.3 

45.7 

15 

48.2 

48.2 

48.4 

50.9 

16 

49.7 

49.7 

49,9 

51.1 

17 

45.9 

45,9 

46.1 

46.2 

24 

69.8 

62.9 

70.5 

63.4 

25 

78.1 

69.7 

79.0 

70.3 

26 

79.1 

70.0 

79.9 

70.6 

27 

72.1 

63,4 

72.8 

63.9 

28 

34 

90.8 

77.7 

91.8 

78.5 

35 

100.8 

86.3 

101.8 

87.3 

36 

102.4 

86.8 

102.3 

87.9 

37 

91.7 

78.6 

92.6 

79.6 

38 

44 

99.5 

86.8 

100.1 

87,4 

45 

110.0 

96.4 

113.0 

98.6 

46 

110.0 

96,8 

111.5 

98.0 

47 

99.7 

88,0 

101.0 

88.9 

48 

54 

97.1 

88.5 

97,5 

88.5 

55 

107,7 

98,5 

108,0 

98.7 

56 

108.1 

99,2 

108.8 

100,0 

57 

98.2 

90.5 

99.0 

91,2 
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TABLE  III 
(continued) 


go  NONLINEAR 

e 

NO  CONTROL       CONTROL 


LINEAR 
NO  CONTROL      CONTROL 


64 

88.7 

83  ,'8 

89.3 

84.3 

65 

99.3 

94,1 

99.5 

94.8 

66 

100.4 

95,4 

101.1 

96.5 

67 

91,9 

87.6 

91.9 

88,0 

68 

69 

74 

79.0 

76.5 

79,4 

77.0 

75 

89.5 

86,9 

89.9 

87.3 

76 

91.3 

88.8 

91.3 

88.6 

77 

84.4 

82,2 

84.8 

82.6 

78 

79 

84 

69.9 

68.8 

70.1 

69.0 

85 

80.2 

78.9 

80.1 

79.3 

86 

82.6 

81.4 

83.0 

81,7 

87 

77.4 

76.1 

77.0 

76.5 

94 

61.9 

61.4 

62.0 

61.5 

95 

71.9 

71.5 

72.1 

71.6 

96 

75.0 

74.5 

75.2 

74.7 

97 

70.9 

70.5 

70.9 

70.6 
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TABLE  IV 


SINUSOIDAL  LINEAR  MODEL  OUTPUT  WITH/WITHOUT  FEEDBACK  CONTROL 


AT  ld  =;45,0  Rad/Sec 
e 


TIME 

CONTROL 

CONTROL  EFFORT 

NO  CONTROL 

T 

Y 

3 

u 

Y 

0.780 

54,38 

4,28 

61.02 

0.785 

34.33 

2.64 

38.03 

0.790 

12,50 

0.86 

13.08 

0.795 

-9.96 

-0,95 

-12.54 

0.800 

-31,91 

-2,7  2 

-37.52 

0.805 

-52.23 

-4,35 

-60.56 

0.810 

-69,85 

-5.76 

-80.49 

0.815 

-83.86 

-6.87 

-96.26 

0.820 

-93.56 

-7.62 

-107.08 

0.825 

-98 .44 

-7,99 

-112,38 

0.830 

-98.23 

-7.94 

-111.88 

0.835 

-92.97 

-7.48 

-105.62 

0.840 

-82.90 

-6.64 

-93.92 

0.845 

-68.60 

-5,45 

-77.37 

0.850 

-50.74 

-3.97 

-56.85 

0.855 

-30,27 

-2.32 

-33.39 

0.860 

-8.34 

-0,53 

-8.21 

0.865 

14.22 

1.29 

17.40 

0.870 

35.95 

3.04 

42.10 

0.875 

55.82 

4.63 

64.64 

0,880 

72.82 

5.99 

83.85 

0.885 

86,06 

7.03 

98.73 

0.890 

94.87 

7.72 

108.53 

0.895 

98.79 

8.00 

112.73 

0,900 

97.61 

7,88 

111.11 

0,905 

91,42 

7.34 

103.79 

0,910 

80,50 

6.43 

91.11 

0.915 

65.44 

5.19 

73.74 

0.920 

47.01 

3,68 

52.56 

0.925 

26.49 

1.98 

28.68 

0.930 

3,94 

0.17 

3.31 
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1,   The  Characteristic  Equation  For  Feedback  Control 
Applied 

The  signal  flow  graph  of  Fig.  2  was  modified  to  show 

feedback  gains  and  is  shown  in  Fig,  11.   Using  Masons 

gain  rule,  the  characteristic  equation  was  derived  from  the 

S.F.G,  (Eq.  4,B.10).   The  resulting  equation  was  compared 

to  characteristic  equation  derived  from  the  modified  system 

matrix  with  feedback  of  Appendix  D.   It  was  noted  that  they 

differ  by  the  terms 


a21a33  and  ~a23a31  ' 

These  terms  very  nearly  cancel  and  are  insignificant  Ca„j 

a„„-a9qa-  = . 5701)  when  compared  with  the  other  larger  terms 
in  the  characteristic  equation.  Therefore,  their  omission 
from  the  signal  flow  graph  of  Fig.  11  is  valid.  The  deri- 
vation follows : 

L        =    —1  = = C4.B.2) 

Is  s             s 

(k0)(-k~)(-k7)  ~kqkfik7                 ~a33 

T        -  3  6  7  _    3    6    7  _    6J_         (4.B>3) 

L2    "                      s  s                              s 

(-k0)(k„)(-ku)(-kft)  -k  k  k  k               -a 

T        -    2         3            4 8_  .  ,  ,2    3    H    8t        _          21          ( 4  ,  B . 4) 

ijq    "                            2  2                               2 

d                            s  s                            s 


L4 


S  s 


(ky)Ck   )(-k4)(^gl)  k^fy  a23 


s 


s  s 
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A  *  1-.[L1  +  L2  +  L3  L4  +  L5]  +  L±    L3  C4.B.7) 


s     s        s  s 

(4.B.8) 


A  =  i  Cs3+g3s2  *  a33s2  +  a21s  +  a23g2s  -  a23§1 


+  a2i§3)  (4.B.9) 


s3+(a33+g3)s2  +  (a23g2+a21)s  +  (a21g3  -  a^g^ 

(4.B.10) 

2 ,   Root  Locus 

The  following  root  locus  plots  are  a  study  of  the 
effects  of  changing  the  three  different  control  feedback  gains. 
Many  weighting  combinations  were  tried  and  it  was  discovered 
that  the  system  was  most  sensitive  to  changes  in  air  flow  rates 
This  observation  was  in  agreement  with  the  analysis  of  Section 
II  that  the  system  was  most  sensitive  to  kR .   Each  of  the 
following  root  locus  plots  and  associated  table  of  values  was 
chosen  as  representative  of  the  alternatives  tried.   The 
characteristic  equation  used  is  shown  as  Eq,  (U.B.IO).   The 
Naval  Postgraduate  School  program  for  the  IBM  360/67  computer 
was  used  to  generate  the  tables  and  plots.   The  program 
listing  is  shown  in  Appendix  C, 
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In  the  first  study,  the  constant  gains  used  were 

SU   =   425'0   >   g12   =   91 

and  s13  wa,s   vaTied  over  three  decades  starting  from  g  -  =  .  1 . 
The  resultant  plot  and  table  of  values  are  shown  in  Fig.  12 
and  Table  IV  respectively.   For  ease  of  reference,  the  sys- 
tem poles  were  given  letter  designations  for  the  characteris- 
tic equation 

(S+A)(S+B)(S+C)  =  0  (4.B.11) 

where 

A  =  (  +  32  .  25-j  35.04) 
B  =  (  +  32 . 25  + j  35 . 04) 

C  =  (2.5xlcT4+jO) 

As  g13  nears  -20,  pole  a  rapidly  enters  the  real  axis  and 
migrates  towards  negative  infinity  and  pole  B  moves  upward 
across  the  real  axis,   Pole  C  moves  in  a  negative  direction 
until  g  _=-26  and  then  moves  below  the  real  axis  tracking  the 
mirror  image  path  of  pole  B.   It  was  seen  from  this  plot 
that  as  g..  q  nears  -20  that  the  complex  pair  undergoes  the 
greatest  movement.   In  this  study  there  was  no  point  at 
which  instability  occurred,   When  the  gains  were  tested  on 
the  system,  there  was  no  noticable  decrease  in  vertical 
acceleration. 

The  next  phase  of  this  investigation  uses  constant 
gains  of 
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TABLE  IV 

Roots  For  Modified  Pressure  Lift  Only  Model 3  g  ^425.0,  g   =91, 

g13  Varies 

k6  A               v  B  ,  C 

-0,099  -31.49+J46.96  -31 .  49-.j  46 ,  96  -1.49+jO 

-0,127  -31,52+j46.88  -31 . 52-j 46 . 8 8  -1.586+jO 

-0.162  -31.52+j46,86  -31 . 52-j 46 . 8 6  -1.611+jO 

-0.207  -31.53  +  J46.84  ,-31 ,  53-j  46  ,  84  -1.643  +  jO 

-0.263  -31.53+J46.81  -31 . 5 3-j 46 . 81  -1.684+jO 

-0,335  -31.54+j46,76  -31 , 54^ j 46 . 76  -1.736+jO 

-0.427  -31.55+J46.71  -31 . 55-j 46 , 71  -1,802+jO 

-0.544  -31.57+J46.64  -31 t 57- j 46 , 64  -1.887+jO 

-0,693  -31,59+j46,56  -31 . 59-j 46 . 56  -1.995+jO 

-0,883  -31.61+J46.45  -31 . 61-j 46 , 45  -2.134+jO 

-1.125  -31.64+J46.30  -31 , 64- j 46 , 30  -2.312+jO 

-1.433  -31,68+j46.12  -31 . 68-j 46 . 12  -2.540+jO 

-1.825  -31.72_j45,88  -31 , 72-j 45 . 88  -2,835+jO 

-2,326  -31.77+J45. 58  -31 , 77- j 45 . 58  -3,216+j0 

-2.962  -31.84+J45.18  -31 . 84-j 45 . 18  -3.711+jO 

-3.773  -31.91_j44.67  -31 , 91- j 44 . 67  -4.359+jQ 

-4.807  -31,99+j43.99  -31 , 99- j 43 . 99  -5.214+jO 

-6,124  -32.06+j43,08  -32 . 06- j 43 , 08  -6,359+jQ 

-7.801  -32.09+j41,87  -3 2 . 09- j 41 . 87  -7.923+jO 

-9.937  -32.03+J40, 20  -3 2 . 03- j 40 . 20  -10,134+jO 

-12.658  -31.70+J37.87  -31 . 70- j 37 , 87  -13.450+jO 

-16.125  -30.60+J34, 50  -30 . 60- j 34 , 50  -19.032+jO 

-20.540  -26 . 95+ j  30.22  -26 . 95- j 30 . 22  -30.645  +  jO 

-26,165  -48.97+jO  -20 . 53+ j 29 . 31  -20 . 534- j 2 9 , 31 

-33,331  -64.98+jO  -16 . 02+ j 3 0 , 86  -16 . 023- j 30 . 86 

-42.458  -80.26+jO  -12 . 84+ j 3 2 . 65  -12 , 83 8- j 32 t 65 

-54.085  -96,55+jO  -10 , 36+ j 34 . 38  -10 . 3 67- j 34 , 38 

-68.90  -115,02+jO  -8,36+36.02  -8 . 363-j 36 , 02 
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g-n  s  -195.32,    g._  =  34.9 


12 


while  g  ,  varied  over  two  decades  from  -.1.   These  gains 
were  calculated  from  the  Riccati,  with  weighting  matrix  j 


Q  = 


and 


10  0  0 
0  122  0 
0       0     .5 


(4.B.12) 


R  =  ,03 


(4.B.13) 


The  resulting  plot  and  table  of  pole  values  are  shown  in 
Fig.  13  and  Table  V  respectively.   As  g"   varies,  the  poles 
follow  similar  paths   as  in  the  previous  study.   The  greatest 
degree  of  pole  movement,  or  sensitivity,  occurs  at  g   =-16.0. 
This  set  of  gains  produced  instability  when  g  _>-1.08. 
In  the  final  study,  the  constant  gains  were 


gll    =  -1688.095,  g12  =  119.16,  g13  =  -18.3 


06 


Again,  g_ ~  was  made  to  vary  over  three  decades  starting  with 
g    =  -,1.   The  resultant  plot  and  table  of  values  are  shown 

-L  o 

in  Fig.  14  and  Table  VI  respectively.   These  gains  were  cal- 
culated using  the  Q  and  R  values  shown  in  (4,A.5).   Again, 
the  poles  follow  similar  paths  but  pass  closer  to  the  origin 
than  in  the  previous  two  studies,  causing  a  greater  net 
decrease  in  the  systems  oscillatory  tendencies ,   As  in  the 
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ROOT    LOCUS    PLOT    FOR    MODIFIED 
PRESSURE    ONLY    LINEAR    MODEL 


Figure  13 
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TABLE  V 


Roots  For  Modified  Pressure  Lift  Only  Model,  g   =-195.32, 

g12=3U?98,  g13=varies  11 

kg  A  B  C 

-0.099  -32,65+j40,60  -32 , 56- j 40 , 60  0,81+jO 

-0.127  -32,66+j40,50  -32 , 66- j 40 , 50  0.71+jO 

-0,162  -32.67+j40,48  ^3 2 . 67-j 40 , 48  0.68+jO 

-0,207  -32,67+j40.45  -32 , 69-j 40 , 45  0.64+jO 

-0,263  -32.67+J40.41  -32 . 67-j 40 , 41  0.60+jO 

-0.335  -32.68+J40.35  -32 . 68-j 40 . 3 5  0.54+jO 

-0.427  -32,69+j40.29  -32 . 69-j 40 . 29  0,46+j0 

-0.544  -32.70+J40.20  -3 2 . 70- j 40 . 20  0.37+jO 

-0.693  -32,71+j40.09  -3 2 . 71-j 40 . 10  0.25+jO 

-0.883  -32 . 72  + j  39.95  -32 . 72- j 3 9. 95  - . 09  + j  0 

-1.125  -32.74+j39,78  -32 , 74-j 39 , 78  -Q.ll+jO 

-1.433  -32 . 76  + j  39.54  -32 , 76- j 3 9 . 54  -0.37  +  jO 

-1.826  -32 . 78  +  j  39.25  -32 . 7 8-j 39 . 24  -0.71+jO 

-2.326  -32.81+J38 .86  -32 , 81- j 3 8 . 86  -1.14+jO 

-2,962  -32 . 84  + j  38.35  -32 . 84- j 38 . 3 5  -1.71+jO 

-3.774  -32  .  86  + j  37.68  -3 2 . 86- j 37 . 68  -2.46  +  jO 

-4.807  -32 . 86  +  j  36.79  -32  .  86- j 36  .  7 9  -3.47  +  jO 

-6.124  -32.81+J35.59  -32 . 81- j 35 , 59  -4.85+jO 

-7.801  -32 , 66  + j  33,93  -32 , 66- j 3 3 . 93  -6.80  +  jO 

-9,937  -32.23+J31.54  -3 2 . 23-j 31 . 54  -9.74+jO 

-12.66  -31.01+J27.93  -31 . 01- j 27 . 93  -11.48+jO 

-16.12  -26.18  +  J22.73  -26 . 18- j 22  .  73  -27  .  87  + j  0 

-20.54  -47.81+jO  -18 . 37-j 23 . 83  -18 . 3 7- j 23 . 83 

-26.16  -61.49+jO  -14?27-j26.51  -14 . 27- j 26 . 51 

-33,33  -74.16+jO  -11 , 43-j 2 8 . 89  -11 . 43- j 28 , 88 

-42,46  -87,51+jO  -92 , 12-j 31 , 0 5  -9.21-j31.05 

-54.09  -102,51+jO  -73,88-j33 .04  -7,39-j33,04 

-68.90  -120.00+jO  -58.68-J34, 90  -5.87-j34,90 
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Figure  14 
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TABLE  VI 

Roots  For  Modified  Pressure  Lift  Only  Model,  g   =-1688,  g   =119 

g13=-18,3(v3ries)  ll  1Z 

kg  A  B  C 

-0.099  -34 . 66  + j  52,35  -34 , 66-j 52 , 3 5  4.817+jO.O 

-0.127  -34.68+J52.29  ~34 . 68-j 52 . 29  4,751+j0.0 

-0.162  -34,69+j52,27  -34 . 6 9-j 52 , 27  4.733+jO.O 

-0.207  -34.70+J52.25  -34 , 7 0-j 52 , 2 5  4,709+j0.0 

-0.263  -34,72+j52,22  -34 , 7 2- j 52 . 2 2  4.680+jO.O 

-0.335  -34.73+J52.19  -34 , 7 3-j 52 . 19  4.642+jO.O 

-0.427  -34.75+J52.15  -34 . 7 5- j 52 , 15  4.594+jO.O 

-0.544  -34.78+J52.09  -34 . 78- j 5 2 , 01  4,533+j0.0 

-0.693  -34.81+J52.02  -34 . 81- j 52 . 02  4,454+j0.0 

-0,883  -34.85+J51.93  -34 . 85-j 51 , 93  4.355+jO.O 

-1.125  -34.91+J51.81  -34 . 91- j 51 . 81  4.227+jO.O 

-1.433  -34.98+J51.66  -34 . 98- j 51 . 66  4,064+jO.O 

-1.827  -35,06+j51.46  -3 5 . 06- j 51 . 46  3.854+jO.O 

-2.326  -35.17+j51,21  -35 , 17- j 51 . 21  3.585+jO.O 

-2.962  -35.31+J50.89  -3 5 . 31-j 50 . 89  3,240+j0.0 

-3.773  -35.48+J50.46  -3 5 . 49- j 50 . 46  2,794+j0.0 

-4.807  -35.70+J49.91  -35 . 7 0- j 49 . 91  2.216+jO.O 

-6.124  -35.97+J49.19  -3 5 . 97-j 49 , 19  1.463+jO.O 

-7.801  -36.29+J48.22  -36 . 29- j 48 . 22  0.472+jO.O 

-9.937  -36.67  +  J46.91  -36 . 67  - j 46 . 91  -0.849  +  jO.O 

-12.66  -37.10+J45.11  -3 7 . 10-j 45 . 11  -2.648+jO.O 

-16.12  -37.52+J42.55  -37 . 52- j 42 . 55  -5.191+jO.O 

-20.54  -37.74+j38,72  -3 7 . 74-j 38 . 7 2  -9.056+jO.O 

-26.16  -36.96+J32.14  -36 , 96- j 32 , 41  -16.11+jO.O 

-33,33  -42,45+j0,0  -27 . 29+ j 23 , 37  -27 . 29- j 23 , 37 

-42.46  -69 , 58+ jO , 0  -18 . 18+ j 27 . 34  -18 . 18- j 27 . 34 

-54,09  -89,56+j0,0  ^13 . 86+ j 3 0 . 54  -13 , 86- j 3 0 . 54 

-68.90  -109.98+jO.O  -10 . 88+ j 3 3 . 10  -10 . 88- j 33 . 10 
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previous  study,  some  values  of  g  „  caused  instability. 
Particularly,  for  g   >-9,9,  pole  C  was  in  the  right  half 
plane,   The  greatest  sensitivity  was  experienced  for  values 
of  g-.  3  near  g  -  =  -30.   Note  that  the  value  of  g  _  chosen 
via  Riccati  will  cause  stable  operation  and  is  also  removed 
from  the  value  of  greatest  sensitivity. 
3 .   Routh  Criterion 

The  Routh  criterion  was  applied  in  order  to  confirm 
the  range  of  feedback  gains  for  stable  operation  obtained 
from  the  root  locus  study.   The  Routh  array  is  set  up  below. 
The  C.E,,  given  in  Eq,  4  .  B  ,  1 7  ,  is  used  with  the  understanding 
that  all  coefficients  are  absolute  values. 


a21  +  a23  g2 


a33  +  g3  (-a23  gl  +  a21  g3) 


^oi   +  ao:  §19  + 


a23gll""a21g13 


21    -23  *12       a33+g3 


-aOQg1  +  a0,g 


23sl     2163 


The  conditions  for  stability  are 


1)   g3  >   -a33 


since  a33  -  64.493,  then 

g   >  -64,493     for  stable  operations 
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(4.B.14) 


2)   a21g3  >  a23gl 


1    a23    3 


(1.B.15) 


„  2267,782 
gl  <  ~11V2U"8'  '  g3 


g2  <  201,616  g3 

g1  <  (201.616K-64.493) 


g1  <  -13002,821 


3)   a2]_  +  a23  g2  >  .( 


a23gl"a21g3 


a23+g3 


-) 


a23gl-a21g3, 
g2  a23 


(a23gl~a21g3)  (a2la33+a21g3) 


g2  > 


a23+g3 


a33+g3 


'23 


"a23gl+a21g3"a21a33"a21g3 
aQQ+g 


g2  > 


33  &3 


'23 


§9  > 


"a3  3gl"*21*3  3 
*33+g3 


'23 


>  - 


(a33gl+a2la33) 
a23Ca33+g3J 


C4.B.16) 
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If  g3  >  -64  the  g2  >    -201.72, 

The  boundary  limits   of  g? ,  and  therefore  g  ,  are 
very  sensitive  to  g_ , 

For  the  specific  control  chosen  in  Section  IVA?  the 
Routh  array  was  used  to  determine  the  stability  limit  for 
the  gains , 
From  Eq(4,B.14) 


*13  >  ^33 


therefore     18 . 3  >  -64,493 
From  Eq(4.B,15) 

*11  >  a^  «13 

therefore     1688  <  3689.581 
From  Eq(4.B.16) 

(a33gll+a21a33} 
g2  >      a23(a33+g13) 

therefore    119  >  -273.95 

It  was  seen  from  this  study  that  the  gains  were  well 
within  the  limits  of  stability.   Lower  gains  were  realitively 
more  stable. 
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V.   CONCLUSION 

It  has  been  shown  that  the  linear  regulator  design  can 
be  applied  to  the  problem  of  vertical  plane  motion  for  the 
captured  air  bubble  surface  effect  ship.   The  linear  pressure 
lift  only  model,  developed  from  the  simplified  equations  of 
motion,  was  shown  to  be  a  valid  model  which  could  be  used  for 
control  analysis  and  design,   An  optimal  control  was  found, 
based  upon  a  quadratic  cost  function,  using  empirically  de- 
rived weighting  factors,   The  pressure  lift  only  linear 
system  showed  a  reduction  in  center  of  gravity  acceleration 

for  17  <oo  <100  radians/sec,   A  C.G,  acceleration  reduction 
e 

was  observed  in  both  the  linear  and   nonlinear  systems.   The 
shape  of  the  frequency  response  magnitude  curves  were  identi- 
cal.  A  13%  maximum  reduction  was  observed  at  to  =45  rad/sec. 

e 

where  peak  acceleration  occurs . 

The  control  gains  derived,  as  shown  in  Section  IV,  are 
within  the  airflow  rates  currently  available  for  the  XR-3. 
The  Riccati  equation  computer  program  on  Ref.  10  proved  to 
be  a  most  useful  tool,  which  facilitated  experimentation  of 
a  large  number  of  weighting  coefficients,  making  an  optimal 
control  design  more  feasible, 

Finally5  the  analysis  shows  that  the  system  is  stable 
over  the  range  for  the  control  gains  selected.   It  was  seen 
from  the  root  locus  study  that  the  limits  of  stability  narrowed, 
that  the  sensitivity  to  air  flow  rates  increased  and  C.G. 
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acceleration  decreased  with  higher  control  gains,   It  was 

seen  that  for  these  higher  gains  the  C,G,  acceleration  was 

not  reduced  and  the  control,  3  ,  increased  for  w  <17  rad/sec. 

u '  e 

This  was  due  to  the  fact  that  the  buoyancy  effects  are  pre- 
dominant in  that  frequency  range  and  the  buoyancy  terms  are 
at  a  constant  steady  state  value  in  this  model. 

Though  an  optimal  control  was  found  based  on  a  particu- 
lar quadratic  cost  function,  many  other  possibilities  exist. 
Different  defined  weighting  factors  can  be  experimented  with, 
and  some  may  result  in  improved  C,G,  reduction.   The  limiting 
factor  in  such  a  design  is  equipment  limitation.   The  choice 
of  weighting  factors,  and  ultimately,  the  control,  are  even- 
tually determined  by  that  fact. 
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VI.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  following  list  of  recommendations  is  presented  for 
consideration  and  further  study, 

A,  A  more  realistic  system  presentation  could  be  accom- 
plished if  pitch  coupling  was  combined  into  the 
linear  regulator  study.   This  would  facilitate  the 
speed  analysis  study  for  the  avoidance  of  wet  decks 
and  undesirable  pitch  accelerations. 

B.  The  passive  control  technique  discussed  by  Boggio 

in  Ref .  6  should  be  examined  further.   The  introduc- 
tion of  a  flexible  membrane  into  the  plenum  chamber 
may  reduce  vertical  accelerations.   This  technique 
has  been  examined  in  model  towing  tank  studies. 

C.  Since  the  mass  of  plenum  air  X3  cannot  be  directly 
measured,  an  observer  might  be  used  in  conjunction 
with  state  X3  or  perhaps  remodel  using  low  order 
model  analysis.   Perhaps  consideration  should  be 
given  to  a  design  using  only  XI  and  X2  as  states  to 
be  used  for  feedback, 

D,  The  linear  regulator  technique  design  approach  should 
be  compared  in  greater  detail  to  the  filter  design 
technique  proposed  by  Rohr  Marine  Inc ,  in  Ref.  3,   A 
fixed  air  flow  rate,  based  upon  a  filtered  observed 
sea  state,  could  be  a  more  realizable  design  in  terms 
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of  the  physical  equipment  limitations*   The  study 
of  fixed  air  flow,  of  Appendix  E,  shows  that  the 
filter  design  approach  may  offer  greater  control 
over  the  entire  encounter  frequency  range,  (see 
Fig,  E-l) 
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APPENDIX  A 

The  simplified  non  linear  system  of  Ref .  1  is  shown  in 
Fig,  Al ,   This  system,  with  the  buoyancy  term  removed,  served 
as  the  basis  for  study  in  this  thesis,   The  linearization 
and  resultant  linear  system  matrix  is  shown  for  the  simpli- 
fied non  linear  system,   The  pressure  lift  only  model  lineari- 
zation is  also  shown. 

LINEARIZATION  OF  THE  SIMPLIFIED  MODEL 

The  equations  of  motion  are  linearized  using  Taylor 
series  expansion  around  the  steady  state  operating  point. 
Heave  acceleration: 
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3Vm   =    -[nk     +   -JLi  cU^S~)a''2J   Yp    (o) 
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The  state  vector  is 
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—  m 

The  output  variable  is 


y    =    3Pb 
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Writing  the  above  equations  in  state  equation  form 


or 


h    -    X2 


2A  p  g      YPh(0)A.2 

x2  =  -cc-Vi-)  ♦  Cfl-^nrr-)]  x. 


yPb(0)Ab 

+  ["(MV~T0l )]  X3 

m 


X3  =  C-YPb(0)[nV^-l(2_^__)l/2]  \^    h 

b 


+  [-YP  <0)[nk  +-2_1( ^_)1/2]_1   ]  x. 

b       q   Pa   2P.(0)      VmT¥)    - 


Xl  =  X2 


X2  =  a21Xl  +  a23X3 


X3  =  a31Xl  +  a33X3 


The  characteristic  equation  was  found  from  the  system  matrix, 
DET(sI-A)  =  0 
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to  be 

s3  +  62s  +  2183s  +  789  =  0 

which  factors  into  j 

(s+.37)Cs+30.82+j34,8)(s+30.82-i34,8) 

TAYLOR  SERIES  EXPANSION  OE  PRESSURE  L1ET  ONLY  VERTICAL 
EQUATIONS  OE  MOTION 

Defining 

W 
Wm  =  M 


F    =  -2. 
17  pm    M 


then  incremental  acceleration  due  to  pressure  i: 


z  (o)'+  az-«g-  cf  co)  ♦  sFpm) 


3Z  =  -  aF^m 

pm 


Defining  the  air  mass  volume: 
a 


V(0)  +  3Vm  =  nQ    -  nk  CPb(0)  +  3Pb) 

m        m     -i'v     ^ 
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Then  the  incremental  volumetric  change  in  air  mass  is; 


»         q    b     p^    2p^(Q)        b 

D 


and  the  incremental  force  due  to  pressure  is 


3F   =  A,  3PV 
p    b  b 


Incremental  air  flow  into  the  plenum 


3q.   =  -(nk  )  3PT 
nm       q   b 


Incremental  air  flow  out  of  the  plenum 


3q     =  -  JSJ:  ( *_)1/2  3PT 

Hout       pa    2p-(Q)        b 


Incremental   gauge   pressure    (plenum) 


3V  3V 

3Pb    =    YPb(0)    l—$J-   jj-fal 
m  d 


Incremental  plenum  volume 


3Vb  =  AbLd 
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Incremental  Pressure  (plunum) 


D        D 


Incremental  Draft 


3L,  =  3Z 
a 


Combining 
Acceleration  ( 3Z) 


32  =  "  T  ^ 


A,  3V      -A,  3Z 

m       b 


b  m 


Volumetric  air  mass  ( 3  V  ) 

m 


ca,   pa   , /9      3vm     ay, 

*      pa      2PTC0)  b  Vm(0T       VbrST 

D 


C  A,         p  n  ,  „  3V  A.  3Z 

m  <1      pa      2PT(.0)  b  Vm(0)         VbT0T 

b 
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Plenum   Pressure 

3Pb  ■  YV0)   lTTST  +  ™  3Z] 

m  b 

A  series  of  system  coefficients ,  k  ,  were  defined  for  ease 

'   n? 

of  writing  and  are  used  throughout  the  thesis.   They  are, 

2W  V        WT  L 
X  w 


,         1 

K2     V  CO) 


k0   =  yP^(O) 

o        D 


A, 
k    =  — 


I 


k„   =  nk   + 


n  1         a 


1  V 


^     pa     T  2P~(0) 

D 


.      _     1 

k7    -  rrsr 

m 


k8   =  Ab 


k  .  =  1  constant 
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APPENDIX  B 

A  brief  step  weight  removal  study  was  performed  and 
compared  to  Ref.  1  as  an  additional  check  on  system  validity 
The  IODE  remote  terminal  CP-CMS  program  on  the  IBM/ 360  was 
used  in  this  study.   A  sample  listing  is  shown  in  Appendix 
C. 


7S 


X  1 


a. '-02 


:-.'* 


.K5 
J 


0.; 


B.  750f-C2 


).[  75 


t.2o; 


n.O  e.75C=-02 


)  a ; 


n.2oi 


<~.;50 

♦    C.500 


J  .  <.  5  I 


♦    0.402 


•    C.353 


t    0.304 


Ic. 


•54 


:.^r>o 


Figure    Bl 
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C.3-<. 


6  -  3  C  1 1.  -0  i 


1  .  1 1 0 


CIS! 


I 


i.  e«         ♦ 


>■    >.  >■ 


U.7V- 


<    -1.6' 


-2.c.3  ♦ 


S.il 


-A.  C2  ♦ 


■6.20  « 

T  0.C 


VlLCCITy    VS    TIMS 


0.1.! 


'•-.180 


♦    -2  .S3 


f    -A. 02 


-5.11 


♦     -'i.21 
._  ^ — 

0.  MO 


Figure    B2 
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-<2  .t 


C.C 


I 


)CCE-02 


I .01OE-C  1 


C.I'jC 


-A3. 9  ♦ 


f     -41.3 


55.1  * 


-  6C-1 

T 


-tl.3  ♦ 


-67  ,6 


I 
-7J.8  ♦ 


0.0  5  .OCOE-02  l.OOOE-Ql 

PLcNtM    MS    FLCW    SiTc    VS    TIME 


0.150 


4    -55.1 


♦    -61  .3 


»    -6''. 6 


4    -7?. 9 


—  ♦  — 

n.2  3C 


Figure   B3 
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lj  r , 


s.  oool:-i;^ 


ifc-Ci 


">.l  'O 


C  .  2  C  C 


137. 


Figure    B4 
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3.v<; 


I  -- 


t       H.~/9 


1  -,, 


:/d. 


I 


f    -3  76. 


0.0  5.CCCr-.V 

ACCEL  F  =  \T  ICN    VS    TINE 


L.cni=-ci 


C-15C 


-505, 


»    -*>?3 


0.<.'C 


APPENDIX  C 

The  following  are  the  listings  for  the  various  computer 
programs  used  in  the  preparation  of  this  thesis. 
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//    ■■< 

/  /  L  J  s 

r  ■.  z 

//oJ  . 

I-    .,   V 

S  Y  S  7 

-22*  7 

-  L  3  i  'J 
B     -' A  T 
^.  ) 
C  U  T  F  L 
-2267 

D      M  \   1 

C    :^AT 

ic  i :  g 

0.  ) 

j .  i 

r 


C  L  T  *  I Z  C 

* .  -:  v  : ' : .   r.  * 

L  L  C  "  ;  Y  S  L  i  E  ( 9  f  r  A  T  T 

:•  v  - .  i  l  c 


c-::n"     l 

.-'    "  i""  ;  *  y 

I .  C 

.73  2    C.T 


r  I  < 

r    v  j 

.  1  c  2 

F  IX 
F  IX 


.0 

I  X 

.c 


u)    ,,  \ r"     «    N  A  " 

.. .  j 


>m 


-11.246 
-64.49  3 

1.0 

-11.246 

2.-") 


*  • 

.  5 
jA  ,(-'■:  c;     pro     td  a  n,,  ct  -at     S'-)LiJTICN 

"  2  5  .  )  j  2  ; 


i  j  i : : 

j  •  • 

CO 
E  L  A  N  K 
BLANK 
B  L  A  f  i  K 


C  A  3D 
C  £  ?  C 
C  A  •  ) 


^  •  ».. 
1 

c.c 


r.'l 
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//       x    r    G  •• 
//"•    L  .  '■■  F!J" 

~*  ri      *  -■    r    : 
r  :■■•  r  g  '  :•     .  p  l  ; ' 


rr    * 

£  i  '•  H  T  L  *  T  Y 


CYC      * 


C  f  "• '    ■      I  /  C  A  R  r  A  /  S »  C ,  fJ  •  G  f  H 

C  f.  M  S  T  -=:2o7.7-2 

Cn mST  GP=2  2. 17* 

P  A  *  A  N  K  6  =  c  .  1  -  A  :  2 

rvMAN t: 

LGGW=f AVFIC. ) 

l*=10.»*LCGfc 

£=CvFL>(r  .  f  v.  ) 

V=(  2.*  2.1  '!55"r"*l  /W 

*=c2  .c  319/V 

K=-A*t  SlM>))/y 

\=-A*<  S**2) 

■I  =  (  S  v  '■  2  )  +  (  1  .■  .  ^G  7  *K  6  *  S  )  +  ?  2  6  7  .  7  c  2 

G  =  N/C 

I-  =  G  /  ( 1  +  G  ) 

-  •  =  -  .  A  I  (  C-  ) 

T  v  =  \  T  ,v  t  Q  (  G  ) 

5->=?     AL(M 

'  I}-  =  AIN£C<I-) 

FH  I  =  *7  .."V?  7M?  (  !»',  '?-  ) 

IF(PM.G7.C. )DH!=P  H!-3f. 

MG=CAES (G ) 

V  lQrl  =  c  fl  a C  (  H  ) 

MC3s20«*4LCG10<  *  A  G ) 
NCBh  =  2w.*flLCGl  MV'AGH) 
G  ^  r  A  L  =  L I  M  T  (  : .  Z  ,  1  *>  ) .  1 ,  P  ?  ) 
CiyAG=LINI-(C.C,A6.Ct  !-V  ) 
PHI  18  >  =  Ph  I  +  18  :">. 

c  w-==-iec. 

C  C  ?  L  =  0  . 

c  i  v.  p  |_  = 

-ALL      D«VlG  (   1  ,  1   ,LCGW.  vGa) 

C  'l'_    CF-InG  (  1  ,2,i  ^w    <*< 
T-P.vinpl 

CALL    JMjSMNFLCT) 


*  - 


1.01, 


,yCa.  10BH,Ph  !,PH I ISO 


L    PINT  I  M2.  J»G'lT  =  .^5,CrlS  =  .,n5 


° ' :  ■ "  ". . :  5  .  ■*  ,  >  r 

ID 

=  •-■;/  <<6  =  12. 0 

F  A  ?.  •-'.  N  K  6  =  1 9  .  J 

"  10 

ST  IP 

//C.SY SPRINT 


»•<--'» 


I '1,"fM  ,  *  AG  ,  VC3H  ,  v\GH 


SVSPUT=Af DSN  = 


//pl:ii  . s  v s:n   gc   * 

CAB    S"f    K6  =  6.14*2<; 


c      1 


PL 


PAG  (Ct  )     VS.    L 


l  J       A 


.-  5     K6=12  .C 


BCOc     FLCTi     M^G(Cc)     VS.    LGG    V, 

CAB    S     3    K6=19.C 

BUDS    FLOT,    MAG(C8)    vs.    log    k 


fiS 


V 


// 

Y 

' 

i 

/  /' 

L  . 

r      r-  - 

* 

-"     T       - 

w 

i 

"  I 

*    ; 

r  ,\    — 

'- 

r    | 

l. 

- 

r  r 

-• 

?  L  " 

"~  =  t 

"  I 

'  -  L 

• 

'  K  - 

r  .  t 

T 

• 

)C  S 

♦  z  2  r: 

r  i 

p  .■  : 

:      i 

V 

=  5 

..  • 

•  •■  • 

--     T       • 

\ 

~  '     ' '  !  L  /  ~  "  r  '         xr 

T  =  .  H  1  5  »     f   ■  L  :  =  .  r  0  r> 


f»-cg=i.cc 

G  =  3  2  .  .1  7  A 

rL  =  (2.r  '.'^KH^'M  . 5)  )  /*. 

- 1-  C  =  .  v  C23  7t 

^B=2C.f  . 

N=5. 

C  N  =  .  r> 

F  A  =  2  1 1  6 . 

Ia=60C7, 

N  =  A  /  G 

- '    i!=((-V./N)  +  (FL/V))/(-«P/v) 

*C  =  1. 

u  *  C  =  ^  i     ' 

CaNMC   I-I-KC*  P3=  I  ) 

L L  =  Q/  (  C?v*  S  Q P  T  U  2  .  * P 8  3  !  ) /P  HZ  )  ) 
V4=3S3  . 
G  A  v  =  1 .  4 
IS-2.5 

^L=?^.^hCM(("aT^')/Pn^(l./GV'M} 

Lw  =  (  ?.*3  .141592654*G  )  A-r 
7F;  T.i=(i.J  415<;2£54*!<M  /LI* 

Hl='-.M(CI'    I>'TM))/Th!'-': 
I  V  ,*.  7  I  V  . 

T  ;    =  '  -  ;  T  ,\  (  ^  -  *  7  I  v  :  ) 

L  C  =  Z  S  +  2 

Vis*yN-(AE*LC  )-(Kl*(TN  n 

FB=P^  <Nd/(\.B*«hQ)  )  <*G  *v 

F 6  6=  P r  -  F  / 

■  ■  F F  c  =  ■;  c r.  ( P  p  c  > 

C  )IT=(CMAL  )  M?  ^~(  (2.**FBB) /°Hn) )*(*PPB/F 331 

C  7  N  =  N  *  (  C  1 C  -  K  C  *  P  9  3  ) 

Z2 CCT =- I * E / N ) * F  3  E+ ( W/ v  } - ( F 1/ M  ) 

■v     -  DT=RPC  »  (CIN-CHUT) 

2 :  :  :'  =  T  IS  T  G  fi  L  ( •: . .    Z  3  n  1  T  )  I 

2  =  IN7GF  K-2.C  ,Z3nr) 

N -•="••  TGP  L  (NBC  »    ^P00T) 

7  ::  N  *.,N  iL 

CALL        ^.J"V^(f  FLCT) 

,   IC 

Fi^,"     V.  -■  =  5  1  . 

=  M0 

P  :. "  i   "    1ft    =39. 

m 

P.'..v.'   '     V     =38. 
r  JC 

P  ;  :•  <  v     ^  '  •  =  3  3  . 

PA: AV     ^  "=31- 

r  ^  a 

CTQp 

//C.SV3PFImt    "  C     S  >  S  C 1 7  =  A, »  0  S  N  = 
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OPTIMAL    CGNTROL/KALMAN    FILTER    POOGPAM 
PROBLEM    IOriNTIFTCATICN    -    h  =  AVS    CTNTROL 

The  a  matrix 

0.0  l.OOOJOOuOD  30       ?.••> 

--2.26778200C  C3       0.0  -i. 124800000  01 

-1.3002bC30C  C4       C.C  -6.44S300C0D  Oi 

Tht  e  MATRIX 

0.0  CO  l.OOCOOOOCO   00 

THE    C    MATRIX 

-2.26773200C   03  0.0  -1.12480000D    01 

»**   CONTROL    OPTION    *** 

THE    R    MATRIX 

3.C0000000C-02 

THE    C    MATRIX 

1.0O0CO003C    05  0.0  0.0 

O.C'  l.OOOOJOOOO   03  3.0 

0.0  CO  5.000000000-01 

STEADY    STATE    SCLUT  ICN 

GAINS 

1.68809469C    03  -1.191156:6C    02  1.830578260   01 

«*♦    CONTROL    OPTION    *** 

THF    R    MATRIX 

3.G0000O0OC-O2 

TH1    Q    MATRIX 

1.  JOOGJJOIOC    05  O.'J  1.0 

O.C  l.COCCCOOOO   03  0.0 

0.0  CO  5.CCC00000D-O1 

INITIAL    CCNOITICNS 

0.0  CO  CO 

0.0  O.J  0.0 

O.C  CO  CO 

TIME  =t       2.CCC0QKCC  CI 
GAINS 
0.0  0.0  0.0 

TIM2  =,       1.9CC00002C  01 
GAINS 
1.638J95J6C  03      -1.191155710  12       1.83057772D  01 

TIMc  =,       1.6CC0CUC40  01 
GAINS 
1.6B8J9469C  03      -1.191156J6C  02       1.830578260  01 

TIME  =.       1.7CCGC0C6C  CI 

1.638C9469C  03      -1.191156C6C  02       1.83057326C  01 

TIME  =,       1.6CCCCCCS0  CI 
GllNS 
1.688G9469D  G3      -1.191156060  ^2       1.83D57826C  01 
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THc    INPbT    TATA    IS 

GRAPH    TITLC    AMO    SCALING    DATA 

ROOT    LCV.LS    PLCT    FTP     "CCIFI'C 
PRESSURE    ONLY    LINfcAR    fCOb.L 

XSCALF  YSCALC  IXUF  IY=IGH  IWICE  IHIGH  IG^IO 

20.0  ?0.0  2  6  7  6  1 

ORDFF    CF    Th?    Ch/ =  AC TER  I  <T  IC    ECLA7ICN 
3 

CGNSTAM    PART    GF    The    CCEFFICIF.MJ     IN    CI?CrNOTNG    CPDER 
O.IOCICO    01    C.l  0.2267HC    CA    S.57010C    CO 

VARIA3LF    FART    G»:    Th=    COEFFICIENTS     IN    CESCENCING    CPD2R 
0.0  0.104S70     32    CO  C.C 

INITIAL    VALLF    Ci;    THs     V^FIAELE 

C.500CC5    01 

NUMBER    OF    C=CAJTS    70    9C    SPANNEC 

1  .       «. 

THFi    SYST-M    PGL~S    4R-. 

***     I    f    S    L(U'i<TST)     ***       WARNING  (V>I7l-    FIX")  ZQATR  2     (IHR     =       66) 

REAL    PART 

3.12570F-C3  J.12570E-Q3         -0  .2  5 139=-"»3 

IMAGINARY     PART 

0.47£2lc    T2         -C.47621F    C2  C.C 

THE    SYST=M    ZFRDS    A*S 

REAL    PART 

O.!)  O.J 

IMAGINARY    P«,RT 

0.0  C.C 

ALL     FGINTS    HAV*    Tr  "    S « N 5    CGCFCIMT5S. 
THIS    PLTT    WILL    \G7    3F    CLTFLT. 
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SPECIFICATIONS 

VARIABLES    6    INITIAL    CONCITIONS: 
XI    =    0,0 
X2    =    C.C 
X3   -    0.0 
T    *    C.C 

CONSTANTS: 
A21    =    2267.7fi2CCG 
A23    =    11.248C0C00 
A21    =    2U6.2*5CCG 
A33    =    10.497C00C0 
B21    =    2267.782000 
621    -    2116.265C00 
632    =     l.CCCCCOCOO 
G    -    32.174C0C0C 
A    =    C.1CCCCCCCCCD    00 
WE   =    A5.5CCCGCCC 
K6    a    6.14422C0C0 
Kll   =    -16S6.C947Q0 
K12    *     119.U60CQ0 
K13   =    -16.2C6CC000 

SPECIAL    FUKTICNS: 

LW    *    <2.*3.141592t5*GJ/WE 

F    '    <SIM62.£319/LWn/<62.8319/LW) 

SEA    -    A*SINtWE*T) 

U2    *    K11*>1*K12*X2+K13*X3 

Y    -    -A21*Xl-A22*X2-B2i*F*SEA 

DERIVATIVES: 
C4XI    /C(T    i    =    * 

X2 
D(  X2    /0(T    )    =    * 

-A2i*Xl-A23*X2-621*F*SEA 
D(X3    /C(T    )    =    * 

-A3i*K6*Xl-A32*K6*X3-B3l*K6*F*SEA+U2 

OLTPUTS: 

TITLE:    CAB    SES     INCREASEC    AIR    FLOW    RATE 
TABULATE:     T  Y 

AT    INTERVAL         0.5C00CO0O0OD-02 
NC    FLCTTING 

ENC    CALCULATION    WHEN    T    .GE .    1.50000 
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-  ——SPECIFICATION- 


VARIABLES  £  INITIAL  CONDITIONS; 

XI  =  0.500000 

X2  =  0.0 

X3  =  -44.5300 

T  =  0.0 

SPECIAL  FUNCTIONS: 

Y  =  2030.26*X1_10.09*X3 

DERIVATIVES: 
D(X1  /D(T  )  =  = 

X2 
D(X2  /D(T  )  =  = 

-2182.68*X1-10.76*X3 
D(X3  /D(T  0  =  = 

-12504. 36*X1-62.0*X3 

OUTPUTS : 

TITLE:   SES 

TABULATE:   T      XI       X2       X3       Y 

AT  INTERVAL    0.300000D-01 
PLOT:   X2 

AGAINST:   T      AT  INTERVAL    0.100000D-01 

END  CALCULATION  WHEN  T  , GE .  0.600000 
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APPENDIX  D 

The  characteristic  equation  for  the  pressure  lift  only 
model  under  control  conditions  was  derived  from  the  system 
matrix  with  the  feedback  control  gains  added  to  the  state 
X^.,   The  result  was  compared  with  the  characteristic  equation 
of  Section  IV.B.l. 
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The  system  matrix  is 


A  = 


'21 


'31 


and  with  feedback  gains ; 


23 


'33 


A  = 


21 


1 
0 


a31+gll    +  gi2 


~23 
a33+g13 


The  characteristic  equation  with  gains  becomes 


S^(a33  +  g13)s2-(a23g12  +  a21)s 


~a23a31~a23gll 


+a2la33+a21g13 


It  should  be  noted  that  the  above  equation  is  identical  to 
that  derived  in  Ref,  1  except  that  the  gains  are  shown  in 
(4.B.4) 

a21  =  -2267.782 

a23  =  -11.248 

a31  =  -13002.8 

a33  =  -64.493 
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From  the  Riccati  equation,  it  was  seen  that  the  gains  had  the 
following  sign  values ; 

~S11  '   +g12  '  "Si 3 
Substitution  of  the  appropriate  signs  yields 


s3+(a33+g13)s2+(a23Sl2+a21)s 


+a21a33+a21g13 


~a23a31~a23gll 


Substituting  in  the  system  value  parameters ,  the  character- 
istic equation  becomes: 


s3+(64,49  3+g13)s2+(2  26  7,7  8  2  +  U,248g12)s-11.248g11 


+2267.782g13  +  .5710 
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APPENDIX  E 
FIXED  AIR  FLOW  RATE  STUDY 

A  brief  study  was  conducted  to  examine  the  effects  of  k„ 

6 

on  the  overall  frequency  response  and  compared  with  the  op- 
timal control  design  results.   Earlier  work  indicated  that 
the  system  response  was  sensitive  air  flow  rate.   Therefore 
three  values  of  k  were  tested  on  the  pressure  lift  only 
non  linear  system  with  the  resultant  response  curves  plotted 

Values  for  k„  are : 
6 

kg  =  6,144      ,  Fig.  El 

kg  =  12.0       ,  Fig.  E2 

kg  =  19,0      m  Fig,  E3 

Note  that  the  high  acceleration  midband  frequencies  are 
attenuated  as  k_  increases.   The  response  curve  flattens 
out  and  midband  frequencies  are  reduced  by  20%.   This  result 
would  tend  to  confirm  the  contention,  of  Ref.  3,  that  fil- 
tered fixed  air  flow  rates  method  of  control  may  have 
advantages.   The  time  response  for  sinusoidal  input  remains 
linear  throughout  the  range  of  interest,   The  time  versus 
acceleration  plot  shown  is  a  time  history  for  an  air  flow 

rate  of  0™  =  60  at  an  encounter  frequency  of  oo  =45  rad/sec, 
^10  v  e 

Appendix  C  presents  the  DSL  Bode  program  listing,   Finally, 
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the  system  peak  frequency  magnitude  plot  for  the  conditions 
of  no  control,  state  variable  feedback  control  and  fixed  air 
flow  rate  control  is  shown  in  Fig,  E8  for  QT  =60. 
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